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grability /Condensed Matter point of view, we simplified inner products of spin chain 

^ Bethe states in a classical limit corresponding to long wavelength excitations above the 

ferromagnetic vacuum. As a by-product we solved a new long-range Ising model in the 
thermodynamic limit. 

Keywords: Integrability in gauge/gravity dualities, Three-point functions in M = 4 
SYM, AdS/CFT correspondence 



ARXlV EPRINT: xxxx.xxxx 



Contents 



1 Introduction 2 

1.1 Formulation of the Problem 3 

2 Classical Limit of the Norm of Bethe Wave Functions (B) 5 

2.1 Derivation 6 

3 Inner Product with a Vacuum Descendent (A) 9 

3.1 Exactly Solvable Case 10 

3.2 Second Toy Model. Regulated interaction 11 

3.2.1 Direct Computation 11 

3.2.2 Path Integral Derivation 12 

3.3 Stochastic Anomaly and the Full A 13 

4 Fixing the Anomaly and the Long Range Ising Model 17 

4.1 Comparing with the Exact Example 17 

5 Classical Tunneling and Conclusions 20 

5.1 Main Result 20 

5.2 Conclusions and Speculations 20 

A Simple Products 23 

A.l A Simple Product in B 23 

A. 2 Non-Stochastic Anomalies 24 

B Path Integral Saddle Point Details 25 

C Four Extra Identities 26 

D Perturbation Theory for the Ising Model 27 

D.l Combinatoric Measure from the Partition Function 28 

D.2 Perturbation Theory 29 

E Generalizations 30 



- 1 - 



Figure 1. A classical extended string decays into two protected classical strings through 
what we call a classical tunneling process. 

1 Introduction 

In recent years there has been outstanding progress in the computation of the planar 
two point functions of single trace operators in M = 4 SYM, see [1] for a recent review. 
Next, one would like to compute all three-point functions. Together with the two point 
functions, these are the building blocks for bootstrapping higher point correlators. 

From the string theory side, the problem can be formulated as a calculation of 
a sphere worldsheet path integral for a string moving in the AdS^ x S 5 background 
with three (non-normalizable) vertex operator insertions. That path integral is not 
yet precisely formulated as very little is known about the exact form of the vertex 
operators. Nevertheless, for two interesting limits there are very promising results at 
strong coupling. The first limit is when two of the operators are classical (with a 
large number of fields) and the third operator is a quantum operator (with a small 
number of fields) [2, 3], see also [4, 5]. 1 The second limit concerns the recent advances 
in understanding correlation functions of three classical string states by computing a 
minimal area surface with three punctures [8, 9], see also [10]. 

In this paper we will consider the correlation function of three classical operators 
at weak coupling. More precisely, we will focus on the case where two of the classical 
operators are BPS while the third classical operator is dual to an extended classical 
string. All the operators are made of scalar fields, that is they are dual to strings with 
non-trivial motion in the sphere. We think of this case as a classical tunneling process 

1 This limit can be generalized to higher point functions of two classical operators and several 
quantum operators [6, 7]. 
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Figure 2. Three-point function of SU(2) operators at tree level considered in [11]. This is 
the simplest non-extremal non-trivial configuration. We consider two protected operators 0\ 
and 03 and one non-protected operator 2 . All operators are classical, i.e. made of a large 
number of fields. The nontrivial contraction (see arrows in the figure) involves L' fields of 2 . 
The non-trivial part of the structure constants is given by the universal ratio *4(u)/£>(u), see 
(1.1). This ratio is in a sense an intrinsic property of the non-protected operator 2 . It only 
knows about the (length of the) protected operators through V which only appears in A(u), 
see (1.3). 

where the extended non-BPS classical string decays into two BPS classical strings, (see 
figure 1 for an artistic depiction). 2 3 

In the rest of the introduction we present the problem we solve. 

1.1 Formulation of the Problem 

For simplicity we will mostly consider the setup introduced in [11] where each of the 
three operators can be embedded in an SU (2) subsector of the theory, see figure 2. On 

2 Thc classical tunneling case is conceptually and technically much more involved than the corre- 
lation function of two classical strings and one quantum protected string recently studied in [2-5] . A 
very important difference is that the result for the structure constants C123 is exponential for the case 
of three classical states and is not exponential for the case of two classical operators and one small 
quantum operator. This is true both at weak and strong coupling. In this paper we will compute the 
exponent at weak coupling. Of course, it would be very interesting to compute the pre-exponent as 
well. 

3 By classical tunneling we mean a tunneling process that is controlled by a classical saddle point. 
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the other hand, it should be clear that the techniques we will develop can be applied 
to a much wider set of examples. Some straightforward generalizations are presented 
in appendix E and a more thorough study will be presented elsewhere. 

In the setup depicted in figure 2, the operators 0\ and O3 are classical BPS oper- 
ators while O2 is a classical non-BPS operator. The corresponding structure constant 4 
^123 wm depend on the number of each kind of fundamental fields in each of the single 
trace operators. It will also depend on the precise form of the operator O* which is 
specified by the Bethe roots u = {ui} or equivalently, by the momenta of the magnons 
X in the Z ferromagnetic vacuum, see [11] for notation and Bethe ansatz review. In 
particular, in the limit where all momenta tend to zero Uj — > 00 the structure constant 
reduces to the structure constant of three BPS operators computed in [12]. 

The tree level result can be written as [11] 5 



Co»o , 
123 1 



u) = cr°° x 




B(u) 



(1-1) 



where L and N are the total length and number of excitation of the non-BPS operator 
O2 while V is the number of contractions between O2 and Oy. The quantities A and 
B are non-trivial functions of the Bethe roots. They are the main focus of this paper. 
The latter is related to the normalization of Bethe eigenstates [13] and is given by 



B(u) 
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Ui — Uj i,j 
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«? + 1/4 
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(1.2) 

while the former comes from the overlap of part of the non-protected operator O2 with 
the protected operator Oi, see figure (2). In the spin chain language it comes from an 
overlap of an off-shell Bethe state with a vacuum descendent [11]. We have 



-4(u) 
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i/2 



u a + i/2 
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The sum in (1.3) is over all the ways of partitioning the set of Bethe roots {u} into two 
groups a and a. The number of elements in the partition a is denoted as \a\. 
In the classical limit [14-16], 



N ~ L — > 00 . 



;i.4) 



4 Thc structure constant is well defined once we normalize the two point functions to one. This is 
implicit throughout the paper. 

5 This formula is valid for finite roots Uj quantized according to Bethe ansatz equations. 
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p(u) and Bethe roots for a two-cut solution with L=168 and N=56 p(u) and Bethe roots for a two-cut solution with L=80() and N=100 




Figure 3. Two examples of real configurations of Bethe roots (pink dots in the complex u 
plane) and corresponding densities p{u) (blue dots above the complex u plane). The depicted 
configurations are the so called two cut solutions. A general classical solution can have an 
arbitrary large number of cuts. The results presented in the main text are valid for all such 
configurations. 

the Bethe roots condense into smooth cuts described by a density p(u), see figure 3. 
Hence, in the classical limit, the structure constant will become a functional of this 
density, 

CS({«i» C?*3°[P]- (1.5) 

It turns out that the continuum limit of (2.1) and (3.1) is quite a challenging computa- 
tion, full of remarkable novel structures. The purpose of the current paper is to perform 
this continuum limit. 

The rest of the paper is organized as follows. In section 2 we study B(u). In sec- 
tion 3 we consider ,4(u). To compute this quantity it turns out that we need to solve a 
long-range Ising model. This is done in section 4 and appendix D. In section 5 we put 
together *4(u) and B(u) to compute C^fp], our main result. In that section we also 
present several comments of more or less speculative nature. Generalizations to more 
complicated three-point functions are presented in appendix E and some complemen- 
tary material is contained in appendices A, B and C. 

2 Classical Limit of the Norm of Bethe Wave Functions (B) 

In this section we initiate our study of the classical scaling limit of Bethe ansatz scalar 
products arising in (1.1). There are two such products that we need to analyse: the 
norm of a Bethe eigenstate (related to B(u)) and the inner product between an (off- 
shell) Bethe state and a vacuum descendent (related to A(u)). In this section we 
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Figure 4. In an interval (u — du/2, u + du/2) one has p{u)du roots described by an approxi- 
mately constant density p(u). As du — > the approximation becomes exact. We can use this 
trick to study exactly quantities which are roughly diagonal as described in the main text. A 
very illustrative example is the toy model (2.4). 



consider B(u) given in (2.1), 

B(u) 

where 

9 
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(2.1) 
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(2.2) 



The goal of this section is to compute this quantity in the classical limit (1.4). The 
final result is given in (2.15) below. 

2.1 Derivation 

The determinant in (2.1) is a determinant of a matrix (2.2) which is a sum of two terms. 
Let us first drop the second diagonal term, i.e. let us consider the determinant of 



M (0) 

i,3 



(Ui - Ui 



+ 1 



(2.3) 



In the scaling limit (1.4), the roots are very large so unless the index % is close to the 
index j the elements of the matrix can be set to zero. In other words, the matrix is 
roughly diagonal which means that the elements of the matrix for which \i — j\ is large 
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are suppressed. To see this we split the roots Ui in blocks. These blocks contain a large 
number of roots each but at the same time, within each block, the density p can be 
approximated by a constant, see figure 4. Then inside each block Ui — Uj = (i — j)/p 
and the matrix elements decays as one over a square distance away from the diagonal. 
This property of the matrix also allows us to say that the different blocks almost do 
not interact with each other, that is 

detM<°> = det ? '~\2~L~\ ~ II det r Xm 2 • ( 2 - 4 ) 

^,J {u i -U j Y + l a n blocks l '° ^ ~ ^ + P" 

The range of i,j in the left hand side is 1, . . . , N while in the right hand side, for each 
block, it should in principle be some large number A, the number of elements in that 
block, which is at the same time much smaller than N. As we will see, the precise 
value for A 3> 1 is irrelevant and the result does not depend on the numbers of blocks, 
provided it is large enough, which is a consistency requirement on this calculation. The 
present analysis leads us to the calculation of the following simpler determinant 

2o 2 

det 4 (2.5) 

ij<A (l - j) 2 + p 2 

for some constant p. Since the matrix is translational invariant, it is diagonal in Fourier 
space. The determinant is simply the product of the corresponding eigenvalues. We 
have 



2p 2 

det — ~ exp 

»J<A (l — ]Y + p z 



f 1 ^ 2p 2 

A / dulog X(p, uj) where X(p, uj) = 7 — 

.Jo _ ' n t^ oo n 2 + p 2 

(2.6) 



The sum can be computed exactly, 

27rpcosh [tt(1-2u)p] 

A(p, u) = — — x • (2.7) 

smh(7rp) 

This solves the problem (2.5). The solution to (2.4) is now straightforward: we simply 
multiply over the several blocks 



det M (0) ~ exp 



/ du p{u) I 
Jc Jo 



du p{u) I duo log \(p(u), uj) 



(2.8) 



where C = C\ U • ■ ■ U Ck are all the cuts along which the roots are distributed and the 
density of roots on those cuts is p(u), see figure 3. 

We now move to the full case, the determinant in (2.1). Since we only have to 
incorporate the term proportional to cV, the full matrix is also roughly diagonal, exactly 
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as in our toy model example. Hence the decomposition into blocks of constant density 
which we performed above still applies. Let us work out the details. The new diagonal 
term in (2.2) contains two pieces: a source and a sum. In the classical limit, the source 
term can be dropped since it is very small 



L 



u 2 + 1/4 



1 

z 



(2.9) 



The sum term in the diagonal part of (2.2) is given by 

2 



(.,. _ 



(ui - u k ) 2 + 1 



(2.10) 



For u k far from U{ we can drop the summands in the classical limit (1.4). For u k close 
to u-i we can use u k — Ui = (k — i) / p where p is the density of the distribution of roots 
around the root U; t . That is, p is the constant density of roots in the block containing 
the root v,;. Hence 



V 



2p 2 



(Ui - u k ) 2 + 1 
Altogether, in each block we have 



XI ,,-1 ,:: 



n z + p 1 



2P 2 



— 2np coth(7rp) = a(p) . 



(2.11) 



det . . 

»J<a \ {i - j) z + p x 



+ a(p)S. 



(2.12) 



We see that the matrix is simply shifted by a multiple of the identity matrix. We 
conclude that the full determinant in (2.1) is simply given by a (2.8) with A shifted to 
a. The integral over oj can actually be performed analytically. 6 The final result can be 
written in the following nice form, 



log det M~ du 2 dp, log [2sinh(7r//)] — p(u) log 
Jc Jo 



sinh[7rp(M)] 
irp(u) 



(2.14) 



This concludes the nontrivial part of the computation of B(u) in (2.1). There is only 
one extra simple product multiplying the determinant in (2.1). That simple factor is 



3 We find, up to factors of in (that is signs for det M), 



log det M ~ / du 



ip(u) / 
Jo 



doo log 



2ttp(u) cosh [tt(1 — 2 uj))p(u)] 



smh[irp(u)] 



— 2np(u) coth[7rp(M)] 



/ — f7rp(u)log \2irp{u)(e 2 * p(u) - 1) 
Jc n v L 



+ Li 2 



(2.13) 
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computed in appendix A.l. It simply cancels the second term in (2.14) so that we end 
up with the neat final result 



p{u) 

log B~ J du J d/ilog[2siiih(7r/i)] . (2.15) 

C 

Of course, we checked this result against explicit numerics from (2.1) and found a 
perfect agreement; see figure 7 for similar numerical checks for the more complicated 
A which we will discuss in the next section. 

Let us end with a small comment. Dropping the source term (2.9) in (2.2) was a 
bit dangerous. Indeed, if we set L = in (2.2), we remain with a determinant of a 
matrix where the sum over the element in each row is zero. 7 Hence the matrix has a 
zero eigenvalue and i3(u)| L=0 = 0. The role of the source term (2.9) is to regulate that 
zero eigenvalue. Indeed, a more careful analysis shows the subleading correction to our 
classical result (2.15) is of order log L/L, coming from the small eigenvalues. 

3 Inner Product with a Vacuum Descendent (.4) 

We will now move to the more involved computation of the classical limit of (1.3) 

qUo={u} u a Ga x ' u a £a,u b £a 

Here we introduced an additional convenient parameter r which we call twist. 8 At the 
end of the day we want to set it to ir. The main new ingredient in (3.1) is the sum over 
partitions of the magnons {u} into two groups: a and its complement a, see figure 5. 9 
The number of elements in the partition a is denoted by \a\. We start by considering 
some simplified toy models to gain some experience and only then move to the actual 
object of interest. 



7 By setting L to zero in (2.2) we mean the coefficient of the source term only, i.e. not changing the 
position of the Bethe roots. 

8 It turns out that introducing the twist is quite useful to regulate several expressions. The use of 
very similar twists to regulate several Bethe ansatz related quantities in the AdS / CFT context can be 
found e.g. in [17, 18] and in the references therein. 

9 In this section we will often use (redundantly) blue for quantities related to the partition a and 
red for quantities related to a. All information is contained in a black and white print-out but we 
believe this section's clarity benefits from being read in a computer or in a color printed copy. 
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Figure 5. The second main ingredient in the construction is the function «4(u) in (3.1) 
graphically depicted in this figure. It is given by a sum over all possible ways of partitioning 
the Bethe rapidities {uj} into the partitions a (blue empty circles) and a (filled red circles). 
For example, for the configuration of seven roots displayed in the figure one sums over the 
2' = 128 way of splitting the Bethe roots into two groups. The purpose of the current section 
is to address the computation of »4(u) in the classical limit, i.e. when the number of roots 
becomes very large. 

3.1 Exactly Solvable Case 

As a first example, lets consider the case where V — in A (3.1). That is, we consider 
the classical limit of 



Al , =0 = y e -M TT " a "st* . (3.2) 

aUa={u} agQ , 
a E a 



Curiously enough this case can be solved exactly before even taking the classical limit. 
We find that 

A L , =0 = (l + e lT ) N (3.3) 

is independent of the rapidities {u}\ The easiest way to see this is to plug (3.2) in 
Mathematica for several values of N and simply check that (3.3) holds indeed. Also it 
is not hard to prove this analytically. First note that At/=o is rational with denominator 
degree that is not less then the numerator degree. The potential single poles of Al'=o 
are where two rapidities, say Ui and Uj, coincide. Such pole comes from two type of 
terms in the sum over partitions: either {iii G a,Uj G a} or {ui G d,Uj G a}. For 
any partition with {ui G a,Uj G a}, there is a mirror partition that only differs by 
exchanging and Uj. The two contributions to the residue at = Uj only differ by 

an overall sign (since — - — = — ). Hence they cancel and the residue is zero. The 

rational function is therefore a constant! To read the value of that constant (3.3) we 
can simply send (one by one) all rapidites to infinity. 

Even though this example looks trivial it will play a crucial role below. 
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3.2 Second Toy Model. Regulated interaction 

As a second toy model example we consider the classical limit of 



e rR w (^4V riFir^- (3 - 4) 

aUa={u} fflSa N ^ a £ a 



a S a 



The only difference between ^4 reg and our object of interest A is that the interaction 
in (3.1) is replaced by in (3.4). This might seem like a harmless 



Ua—Ua 

difference in the classical limit where the rapidities Uj are very large. Indeed, for two 
well separated roots, 



U a - U a + i U a - Ua + 2 

Ua Ua U a Ua 7j \U a U& 



exp ( ) for \u a — Ua\ > 1 (3.5) 



2 

Therefore, one may naively think that in the classical limit A Teg and A are the same. 
However, this is not true. The point is that we also need to consider the interaction in 
the UV region where the roots are very large but the difference of rapidities is of order 
1 so that (3.5) can not be used. As we will see below, this UV region is harmless for 
the toy model A reg but needs to be taken into account for A. Basically the reason is 
that the interaction in *4 reg is antisymmetric but it is not so for A. See appendix A. 2 
for some simple illustrative examples of why this makes all the difference. 

We will now compute (3.4) in the classical limit using two different methods: a 
direct one and the path integral one. 

3.2.1 Direct Computation 

The direct computation is straightforward. First we notice that because of the anti- 
symmetry of the interaction, we can let the product over a £ a go over all the roots: 

^^.ra" n ^4- w 

a\Ja={u} aSa v a £ a 

1 < b < M 

Next we notice that this sum over partitions can be combined into a single product 



A re g — J^J 



MY / is. V M . , ' 

Ua ~ 2 \ TJ u a — u b + t/2 



a=l 



^(^d) II 



u a + \) lh U a -U b ~ i/2 



(3.7) 
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This is then straightforward to write in the classical limit by using (3. 5). 10 We find 

A Teg ^ J dup(u) log l+exp(i$(u)) (3.8) 

with 

d(u)=T-— + -fdv-^ (3.9) 
u J u — V 

where the slash in the integral stands for principal part prescription. 
3.2.2 Path Integral Derivation 

Next, we will re-derive (3.8) by a path integral method. We rewrite the sum over 
partitions in (3.4) as a path integral over densities p a and p^ of particles in the partitions 
a and a respectively. They are not independent since their sum is fixed to be p. Hence 
it is useful to define the single independent variable Q as 

Pa = , Pa = ■ (3.10) 

Then sums over partitions are replaced by 

summand — > / \DVt\ integrand (3-H) 

aUo=u 

where 

log fJ,[Q] = J du[p log p - p a log p a - p a log p a \ (3. 12) 

is the usual entropy factor that accounts for the number of microscopic partitions 
corresponding to a given macroscopic densities p a and p a . 

In (3.11), the integrand appearing in the right hand side is the continuum limit 
of the summand arising in the left hand side. The latter depend on the Bethe roots u a 
and ita while the former is a functional of the densities p a and p & . The summand is 

— =n«" M y"nr^ (313) 

aca ^ a ' 2 / r- a a 2 

ufc " a 6 a 



a £ a 



which becomes 



f i \ lL' . . f , Pa(v) 

integrand = exp / du p a (u) h it + % -j dv 

Jc V u Jc 



u — V 



(3.14) 



10 It is very simple to see that there is no contribution from \u a — U},\ = 0(1) because the interaction 
is anti-symmetric, see the next section 3.2.2 for more details. 
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Putting all the pieces together we find the following path integral representation for 
our sum, 

Ae g = J [VQ}e s ^ (3.15) 

where the action 

S[fl] = log/i[fi] + log(integrand) (3.16) 

reads 

S[Q] = J du [ipaj + p log p - p a log p a - p a log p a ] (3.17) 

where the quasi-momenta is given by (3.9) and the densities are related to the single 
independent variable fl(u) through (3.10). We can now compute (3.15) by saddle point. 
The saddle point is valid because the action is of order L and L — > oo. The saddle 
point condition turns out to lead to a very simple equation for Q which can be easily 
solved leading to 

n(u) = Q c {u) = zp(w)tan^-^ . (3.18) 
Plugging this into the action leads to 

logAreg — S[Q C ] = J dup{u)\og 1 + exp(^(w)) (3.19) 

in perfect agreement with our first calculation (3.8). 

It is curious to notice the following feature of the path integral derivation: at the 
saddle point Vt{u) in (3.18) is purely complex so that the densities p a and p a given by 
(3.10) are not real. This is of course a typical feature of a saddle point method: the 
saddle point can in principle be anywhere in the complex plane. 

This concludes our study of the preliminary toy models. A few other interesting 
simplified examples are presented in appendix C. We will now move to the computation 
of the real object of interest, ^4(u). 

3.3 Stochastic Anomaly and the Full A 

To compute A using the path integral method introduced above, we have to compute 

u a -i/2\ yj u a -u b + i 



^^^nfe® n 




where the average is over all microscopic configurations a and a for fixed macroscopic 
densities p a and p a - We separate the interaction term into symmetric and antisymmet- 
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ric parts by rewriting the integrand as 

V 



3 ir\a\ 



n 



u a -i/2 
u a + i/2 



n 

u a £a,u b ^=u a 



U a -U b + l 

u a -u b -i 



1/2 



n 



u a £a,u b £a 



1 + 



(u a - U b f 



1/2 > 



Pa i Pa 

(3.21) 

where the product in the antisymmetric part can be extended to all roots since the self 
interaction part gives one in this case. Note that the first two products can be simplified 
using the standard expansion (see previous section and appendix A. 2 for more details) 



i/2 



u a + i/2 



n 



u b + i 



ll/2 



u a — u b — i 



exp {ii{u a )) 



(3.22) 



where ^ is given in (3.9). 

We notice that all terms in (3.21) except for the last one are not sensitive to 
different microscopic realization of the densities. This means that the weight for two 
configurations which differ only by a local exchange of one root in the set a with one 
root in the set a is the same. 11 This allows us to write the first two products in (3.21) 
through the densities 



fi[Q] integrand = exp 



iPadu 



n 



\u a £a,u b £a 



-,l/2\ 



1 + 



(u a - U b f 



(3.23) 



Pa) Pa 



while the last product needs to considered more carefully. This remaining average 
is very local since it goes to 1 quite fast. Hence, the nontrivial contribution comes 
from the blocks of roots such that Ui — Uj = 0(1). We refer the reader to appendix A. 2 
where we analyze a very similar, but simpler, quantity. What follows mimics closely the 
treatment of section (2.1) where we separated the Bethe roots into blocks of constant 
density, see figure 4. Inside each block we have — Uj — (i — j)/p- Each block has 
A ^> 1 roots. Of those, —A belong to partition a and —A belong to partition a. 
We must sum over which roots are in a and which roots belong to a. That is, the 
contribution of each block is given by exp ( ) where 



2= lim ^ log y ] [ 

{ni=0,l} a,b=0 



"q(l — "b) 
2 



1 + 



[a - by 



(3.24) 



with 



11 Indeed from (3.22), the ratio of the weights for the configurations that differ by exchange u a with 
U a i gives exp (—i$(u a ) + i$(u a ')) ~ 1 + c ^ff^- where c ~ 1 and we use that u a i = u a + <k ^ k and that 
dui - l/L. 
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1/p 



£ 

{n a =0,l} 
52 "a fixed 



<- 

<- 



+ 



+ 



Figure 6. The effective action is a function of the macroscopic densities p a and p^ only. The 
interaction between the two set of roots however depends on the microscopic way the roots 
are distributed into the two groups of fixed macroscopic densities. That dependence result in 
a contribution to the effective action that we call stochastic anomaly. It is computed by an 
Ising model partition function with a long range Coulomb like interaction log(l + p 2 /n 2 ). The 
figure illustrates the summation of all microscopic spin configurations with a fixed macroscopic 
density in the Ising model partition function (3.24). 



A 

Vn i = ^A (3.25) 

tr p 

held fixed. Note that since p = p a + p ( - x we have Z = Z[p a , p&]. We see that to 
compute this contribution we need to find the partition function of a long-range Ising 
model (see figure 6). Spin up/down (n a — l/n a — 0) is associated to the partitions 
a. I a that contain the root u a , see figure 6. The interaction energy is given by f(n) = 
log(l + p 2 /n 2 ). Finally note that the Ising model (3.24) contains automatically the 
entropy contribution (3.12) and reduces to it when f(n) — > 0. 
The full path integral action then reads 

S[n] = J du [Z(p a (u),pa{u)) + i4{u)p a {u)} (3.26) 

where (3.10) and (3.9). To conclude the derivation of the full >4(u) the main miss- 
ing ingredient is the computation of the partition function (3.24). It turns out that 
computing this partition function is a beautiful but challenging problem. In the next 
section we will show that 

Z(p a {u),Pa(u)) = F(p)-F(p a )-F(p a ) , F(p)= I" dp\ogsmh(np). (3.27) 

Jo 

For now let us take this result and conclude the derivation of ^4(u). In fact, at this 
point, the saddle point computation mimics the end of section 3.2.2 closely. We write 
the saddle point equation following from (3.26) and (3.27), solve it for fl, plug it into 
the action and massage a bit the final result. The details are presented in appendix B. 
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Figure 7. To check the main result (3.28) we compute log.4(u) from (3.1) for the single cut 
solution with mode number n = 1 and filling fraction a = 1/10 for several V /L. For each V 
we compute log^l(u) for more and more Bethe roots (black dots in the bottom right corner). 
We then fit those points (dashed orange curve in the bottom right corner) to obtain the 
asymptotic classical value (solid orange curve in the bottom right corner). The asymptotic 
results for different L'/L are represented by the orange large dots in the main plot. The 
small blue dots correspond to the analytic prediction (3.28). To make the plot easier to read 
we represented the analytics and numerics for different (but overlapping) regions. In the 
overlapping region, the match between the analytics and the numerics is perfect. 



At the end of the day we obtain the very elegant expression, 

q(u) 

log A ~ j J d\i log(l + e^) (3.28) 

uc fe o 

where the contour integral is taken around the cuts and q(u) is a sort of quasi-momenta, 
denned in terms of the density by 

q(u)=T-—+ [ dv^-. (3.29) 
u J u — V 

The results (3.28) and (2.15) are the main results of this paper. The fact that (3.28) 
is so beautiful and simple suggests that there should to be a more straightforward 
derivation. 
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Since the derivation of the main result (3.28) is quite non-trivial it is important to 
check the final result numerically. A perfect agreement is found, see figure 7. For more 
details on the single cut solution used for these numerics see [15] and the appendix C 
of [11]. 

4 Fixing the Anomaly and the Long Range Ising Model 

In this section we attack the only missing ingredient in the derivation of the classical 
limit of -4(u), namely the computation of the Ising model partition function (3.24). 
The result was anticipated above and is given in (3.27). 

The derivation is divided into two steps. First, to get inspiration we performed a 
low density expansion in appendix D. The outcome of the perturbative analysis it quite 
remarkable and iluminating. We find that up to 0(p 9 ) we have 

Z{p a , Ps ) = F(p) - F{p a ) - F(p & ) (4.1) 

with 

F(p) = plog^ + ^- ^ + ^ + G(p 9 ) . (4.2) 

KHJ ^ e 18 900 19845 KH 1 y ' 

The structure (4.1) is highly non-trivial and very inspiring. It leads very naturally to 

the conjecture that (4-.1) should hold for finite density. Assuming this, we are left with 

a single function to fix! This function is then fixed by requiring that the path integral 

for the exact example (3.1) gives the correct result as we now explain. 

4.1 Comparing with the Exact Example 

In this section we will fix the anomaly (3.24) which enters in (3.26). Recall that at this 
point we are after a single function of a single variable: the function F(p) in (4.1). 

The trick is to use the exactly solvable example from section 3.1. This toy model 
is related to the full A(u) by setting L' = 0. The result is given in (3.3), 

log^ =0 = iVlog(l + e' T )- (4.3) 

In particular it does not depend on the distribution of Bethe roots at all. On the other 
hand, if we were to derive it using a path integral we would write down (3.26) with a 
very minor simplification, namely V = 0. That is 

S = Jdu(±(p + n){(? + T) + F[p] - F[(p + n)/2] - F[(p - 0)/2]) (4.4) 
c 
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where 

G(u)=(dv^_ Q( u ) = 1 dv-^- . (4.5) 

J u — V J u — V 

c c 

At the saddle point, S c = S[Q C ] should equal (4.3) for any density p such that 

p{u) = N. (4.6) 



c 

We see that, from the point of view of the path integral derivation, the trivial result 
(4.3) is not obvious at all! It turns out that requiring (4.3) fixes the anomaly completely 
as we will now show. 

We can simplify our life by choosing a very large imaginary twist r = it, t — > +oo. 
In this limit (4.3) reduces to 

log^o-iVe"* . (4.7) 

For large imaginary twists the factor e* T ' a ' becomes very small so that the density p a 
will be very small. This suppression leads to major simplifications in the derivation. 
For example, it translates into Q = —p + 2e where e is very small. Hence 

S = J du (-te - F(e) + eF'(p) + ieQ) (4.8) 
uc 

and the saddle point equation is simply 

-t + i$-F'{e)+F'(p) = (4.9) 

From the result (4.2) of the previous section we know that F(p) = p\og(np/e) + 0(p 3 ). 
Using that we deduce from the saddle point equation that 

6 = I e -* 6 W(P) . (4.10) 

7T 

Finally, on the saddle point, the action (4.8) simplifies to 

S c = J edu= e -^- J e^ +F '^du . (4.11) 
uc uc 

This should be simply equal to (4.7). Thus we have to find a function F such that for 
any density p: 

J e^'^du = tiN . (4.12) 

uc 

This very nontrivial constraint fixes F(p) uniquely. 
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To see this we will perform a small p expansion. We use the structure of the 
perturbative expansion of F(p) deduced from our analysis in Appendix D, 

e F ' {p) = c lP + c 2 p 3 + c 3 p 5 + . . . (4.13) 

Actually we know the first few terms of this expansion from (4.2) but we will not use 
them. Instead we will derive all the c& by imposing (4.12). Expanding the integrand in 
the left hand side of (4.12) in powers of p we get 

J (c lP + i Cl Gp + c 2 p 3 -^Y~ +--^)du = irN (4.14) 

uc 

• Note that the first term gives perfectly tcN if we set c\ = ix. 

• All higher orders should give simply zero. To show that the quadratic term is 
zero we observe that AiirGp = G 2 (u + iO) — G 2 (u — iO) and thus 

J iixGpdu = G{ufdu = (4.15) 

where the integration is done by deforming the contour to infinity. 

• The cubic term can be treated similarly by noticing that § G{u) 3 du = and 

G\u + iO) - G 3 {u - iO) _ tt'V vpG 2 
m ~ ~6 2~ ' 

This allows us to write 

= J (c 2P 3 - ^) du = j [c 2 - p 3 du (4.16) 
u c uc 
since this should hold for an arbitrary p we get c 2 = 7r 3 /6. 

One can continue to apply this strategy repeatedly fixing all the coefficients in this way. 
We find 

n 2n-l 

^p^T)!' (417) 

which leads to the main result of this section, 

e F ' ip) = sinh(Trp) . (4.18) 
It is very easy to check that this is indeed the correct result leading to (4.12): 



J e l(? sinh(7rp)^ = i J ( e lG ^ - e lG ^) du = i j e lG du = ttN 



(4.19) 



uc uc uc 

At the last step we use that G ~ N/u for large u. The result (4.18) leads to (3.27) 
which was the only missing step in the derivation of the classical limit (3.28). 
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5 Classical Tunneling and Conclusions 



5.1 Main Result 

We can now put together the results (2.15) and (3.28) to get the structure constant 
C°*3 (1-1) i n the classical limit. This structure constant involves three classical states. 
Two of them are protected (Oi and (9 3 ) and one is non-protected ((^2)- We call this 
process Classical Tunneling; an artistic depiction is represented in figure 1. In the 
simplest setup - represented in figure 2 - all operators are SU (2) operators. For this 
setup, let us present the final result as a ratio between this structure constant CJ^f and 
the known [12] structure constant involving three protected operators with the same 
R-charges C°23 (times some simple combinatorial factors), 12 



This is the main result of the current paper. Let us summarize the meaning of the 
several quantities in this formula. The non-protected operator 2 is parametrized by 
a set of N Bethe rapidities Uj. In the classical limit they condense into cuts see 
figure 3. The density of roots in these cuts is p(u). The number of Wick contractions 
between operator O2 and 0\ is V which is an integer completely fixed by the lengths 
of the three operators, see figure 2. The quasimomenta q(u) is given in terms of the 
density and V as in (3.29). 

5.2 Conclusions and Speculations 

We conclude with a few comments 

• As discussed in [11], the result (5.2) is highly suggestive. It resembles some sort 
of phase space integration. Recall that the number of roots in each cut are the 
action variables of the theory. Our final result is given by a bunch of integrals over 
the several cuts. Could we shortcut the current derivation by working directly 

12 The non-protected operator O2 contains N scalars X and L — N scalars Z. The number of Wick 
contractions between operator 2 and 0\ is L' which is an integer completely fixed by the lengths of 
the three operators, see figure 2. 



r = 




(5.1) 



Combining (2.15) and (3.28) we get a remarkably simple result 



r = exp / dt 




(5.2) 
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in the action/angle variable basis? Also, the integrand in (5.2) have some kind 
of fermionic flavor to it. Is our computation dual to some fermionic partition 
function? 

• We could also speculate about inspired guesses for the strong coupling value of 
the structure constants. The idea is that the objects arising in the weak coupling 
final result (5.2) have a very natural counterpart at strong coupling. For example, 
the quasi-momenta (3.29) should be related to the eigenvalues of the monodromy 
matrix constructed from the sigma model Lax pair [15]. It is therefore very 
natural to conjecture that (5.2) should describe the strong coupling result up to 
some very minor modifications of the sort described in [11, 15]. Furthermore, 
in the so called Frolov-Tseytlin limit [19], the strong and weak coupling quasi- 
momenta are actually the same [15]. Hence, in this limit, one would expect that 
no modifications are needed. It would be very interesting to try to check this 
conjecture. 13 For more speculative discussions on deforming (5.2) to non-zero 
coupling see [11], [4]. 

• Perhaps most interestingly, let us note that the integrals in (5.2) resemble strik- 
ingly the strong coupling expressions of [8, 9] obtained using completely differ- 
ent techniques. Unfortunately, these results are still incomplete rendering the 
comparison very hard. 14 It would be extremely interesting to complete these 
computations and perform a comparison with (5.2). 

• In this paper we studied the quantities A(u) and B(u) in the classical limit. 
These quantities are the fundamental building blocks for studying the three-point 
functions involving one non-protected operator and two protected operators. For 
simplicity, we focused on the SU{2) setup represented in figure 2. It is simple 
to generalize this setup to more general configurations as illustrated in appendix 
E and in [23]. Moreover, recently the structure constant of three non protected 
SU(2) operators at leading order [11] was written in a determinant form [24]. It 
would be interesting to try to use that simplified form for computing the most 
general structure constants of non protected operators, see also [23]. 

13 At a technical level, it is understood why such match is present for the spectrum problem, see 
[15, 21, 22]. At a physical level, there is no good argument why this match should occur and in fact we 
know that it breaks down at three loops for the spectrum problem. Very recently, it was demonstrated 
that a match also occurs for three-point functions of two large operators and one light operator [4, 20]. 
This match might extend to three-point functions of three classical operators like the one considered 
in this paper. 

14 For example, [8] only takes into account half of the full result, the so called AdS part; the sphere 
part is still missing. Also, in [9] the contribution from the vertex operator insertions is not accounted. 



- 21 - 



The classical results we found might be also interesting from a condensed matter 
perspective. In the condensed matter literature the classical limit (1.4) which 
we consider was first proposed in [16]. It describes long wavelength excitations 
around the ferromagnetic vacuum. There are infinitely many configurations in this 
limit, described by any number of cuts as illustrated in figure 3. We managed 
to compute two important quantities in this classical limit for all possible such 
states. The first quantity is the norm of two on-shell Bethe eigenstates, 



(0\C(u 1 )...C(u N )B(u 1 )...B(u N )\0) 



(5.3) 



where the Uj's obey the Bethe equations in the scaling limit (1.4). This quantity is 
proportional to our B(u). The second quantity is the inner product of an off-shell 
Bethe state with a vacuum descendent, 



(0\(S_fB( Ul )...B(u N )\0) 



(5.4) 



where the u^s are unconstrained complex rapidities in the classical limit (1.4). 
This quantity is proportional to our A(u). 15 Inner products like (5.3) and (5.4) 
appear abundantly when studying correlation functions and form factors in In- 
tegrable models. However, normally, exact results for these quantities are very 
bulky. We hope that the computations presented in this paper can motivate the 
study of spin chain correlation functions and form factors in the scaling limit 
(1.4). Could correlation functions in this limit be measurable in ferromagnetic 
one dimensional chains in the lab? 



Finally, let us mention a by-product of our analysis. It turns out that to simplify 
the computation of the inner product of a Bethe state with a vacuum descendent 
one needs to solve a very interesting auxiliary problem: a long range Ising model, 



Z = lim — log 

L— >oo Zj 



with fixed 




^2 Ua log 



a,b=l 



1 L 



1 + 



by 



'1 - n b ) 



a=l 



Pa 
P 



(5.5) 



(5.6) 



15 Of course, the precise relation between £>(u) and -4(u) to (5.3) and (5.4) depends on the choice of 
normalization of the Algebraic Bethe ansatz creation and annihilation operators. For precise normal- 
ization conventions and precise formulae relating B(u) and A(u) to (5.3) and (5.4) see [11]. 
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To our knowledge this model was not considered before. Quite surprisingly we 
found that in the thermodynamic limit where L — > oo the problem is exactly 
solvable 

Z = f( P )-f(Pa)-f(p-Pa) With f( P ) = Ll2( 2 e ^ P) +^- P l0g(2)-^ 

Technically, this is the most involved computation of the current paper. It would 
be very interesting to see where this exact result fits in the landscape of integrable 
models in statistical physics and to explore the physics of this Integrable long 
range spin chain. 
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A Simple Products 
A.l A Simple Product in B 

To compute £>(u) (2.1), we need to add to (2.14) the piece 

Ui — Uj + i 1 v-^ , 



-'extra 



log J] 



Uj — U-i 



1 + 



1 



(Ui -Uj)'* 



(a.i; 



3 ~ i^j 

The computation of this term is pretty much the same as the computation of (2.10). 
We can again split the roots into blocks of constant density, see figure 4. Inside each 
block we can approximate the right hand side by 



A £ log 



1 + 



smh \np\ 

A log so that 

np 



du p(u) log 



sinh[7rp(M)] 
itp{u) 



(A.2) 

We see that (A.2) exactly cancels the second term in (2.14) so that we end up with the 
neat final result (2.15). 
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A. 2 Non-Stochastic Anomalies 



Consider the following two products (the product is only over j) 

N , • N I • In 

= II , V 2{Ui) = M — A.3 

Ui — Uj Ui — Uj — i/2 

where the ttj ~ N ^> 1 are distributed according to some smooth density p{u). In 
this classical limit two regions can contribute: A UV region where Ui,Uj ~ iV but 
\ui — Uj\ = 0(1) and a IR region where the roots are widely separated. The total result 
is the product of the two, 

T^a — {PaJiR (^ ? a)uv ' 

The contribution from the IR region is the same for both and simply gives, see (3.5), 

(^i)m = (^2) m = exp j where v = u t (A.4) 

and the slash stands for principal part integration. In the UV region we can use 
Ui - Uj = (i-j)/p to get 

( ^ 2)uV= (i-j)/p-i/2 = 11 ^J2 = 1 > (A ' 5) 

where A is some large number (but much smaller than N) governing the size of the UV 
region. As we see its precise value is irrelevant. Note that we get no contribution from 
the UV region simply because the interaction is anti-symmetric. On the other hand, 

(^i)w= n ^ = ? i 5 h <M (A . 6) 
- LJ - , n np 

We see that the classical limit of V\ is more subtle from that of Vi- For the latter, 
we can simply consider the IR region. For the former we need to take into account a 
short distance contribution (A. 6) which we call non-stochastic anomaly. Note that to 
compute this anomaly all we need is the symmetric part of the interaction. 16 

To compute A(u) we see that we need to compute a similar but more involved 
quantity given by (3.26). There we need to take into account the UV contribution 
but the roots are not distributed according to smooth distribution as in the previous 

16 That is we could very well replace (ut — Uj + i)/{ui — Uj) by ^/l + — uj) 2 when computing 

the anomaly. We will do precisely this when computing a more complicated anomaly in the main text, 
called stochastic anomaly. 
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example. Instead we have to average over all possible distribution of roots into two 
partitions a and a. Hence we denote this contribution as a stochastic anomaly to 
distinguish it from the simplest version above. 

There is another closely related quantity called anomaly in the literature [27], [17] 
which is roughly 17 the leading 1/N correction to the trivial UV contribution for "P 2 - 
The anomalies discussed in this paper are not finite size corrections. Instead they are 
order 1 effects that can never be discarded. This is a major difference compared to the 
anomalies in [27], [17]. 

B Path Integral Saddle Point Details 

In this section we present some of the details concerning the saddle point computation 
leading to (3.28) from (3.26). Using the result for the anomaly (3.27) in the action 
(3.26) we obtain the path integral action 



where ^ is given in (3.9) and F(p) = Jq dp logsinh(7r / u) comes from the computation 
of the anomaly. The saddle point equation reads 




) 




uc 



H = fog 



sinh[7r(p + Q c )/2] 
sinh[7r(p-^ c )/2] 





Plugging it back into the action, we get 




(B.3) 



uc 




(B.4) 



g(u) 




(B.5) 



uc 



o 



which is (3.28). 



17 with i/2 — > i in the definition of Vi 
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C Four Extra Identities 



In this appendix we discuss a few more simple cases related to the quantity (3.1) 

v 



A(L',r\u)= J2 < JTH II 



u a - i/2 



n 



u a -u b + i 



We have 

off-shell 



on- 


shell: 


on- 


shell: 


on- 


shell: 



.u a + i/2j u a — Uh 

aUa={u} u a &a N ' u a £a,u b £a 



^(L',7r|u) = , forO<L' <iV 
A(L',t\u) = A(L-L',-t\u)* 

A(L', tt|u) = , for < U < N or L - N < U < L 
A(L,t\u) = (1 + e* 



(c.r 



(C.2) 
(C.3) 
(C.4) 
(C.5) 



The first identity, (C.2) is valid for any set of self-conjugate complex Bethe roots. The 
last three identities hold provided the roots Uj obey Bethe equations for length L, 



uj + i/2 
Uj — i/2 



N 

n 



Uj — life + i 



■Uj -u k -z 



In what follows we shall derive/discuss the four identities (C.2)-(C.5). 

We start by their derivation. The first one, (C.2), is the only identity which is 
not trivial to check. The best is to prove it by exhaustion using Mathematica or by 
using some residue arguments as in section 3.1. Next (C.3). For Uj obeying the Bethe 
equations, we have 



A(L',t\u) 



E 

aUa={«} 



,vr\a\ 



n 



u a &a 



u,, 



i/2 



u a + i/2 



L'-L 



n 



u a 6a 



Un 



i/2 



U a + i/2 



n 



U n 



lift + % 



U a - U b 



n 



Ur, 



u b 



U n 



U b + % 



which is clearly the complex conjugate of A(L — L', — r|u). 18 The identity (C.4) is a 
trivial consequence of the previous two identities (C.3) and (C.2) once we note that 
t = 7T and r = — ir are the same for (C.l). Finally, the last identity (C.5) follows from 
(C.3) and (3.3). 

Let us now move to the discussion of the physical meaning of these identities. 



18 Recall that the set of complex roots is self-conjugate so that {uj} = {uj}* . In particular, in the 
sums, we can replace Uj by u* at will. 
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• As discussed in the conclusions, the quantity A(L',n\u) is proportional to the 
inner product between a vacuum descendent and a Bethe state. These states we 
have iV spins down and V — N spins up. Clearly it does not make sense to have 
L' < N. The first identity (C.2) is a manifestation of this impossibility. Note 
that the inner product interpretation only makes sense when the twist r = n and 
indeed we only get a zero result in (C.2) for this particular value of the twist. 

• The second identity (C.3) resembles strongly the bosonic duality of [17]. Would 
be interesting to explore this analogy further. This might help elucidate the 
meaning of the twist r used throughout this paper. So far we introduced it as a 
mathematical tool. 

• The third identity (C.4) is quite interesting and simple to understand from the 
three-point function point of view. Let us think of the three-point function as 
(0i 1 3 1 (9 2 ). That is we think of the operator 3 as a spin chain operator whose 
average we compute on the states 0\ and 2 as in [25, 26]. Now, 2 is a Bethe 
eigenstate so it is an SU(2) highest weight state of spin (L 2 — N 2 )/2. The state 
0\ is not an highest weight. Instead, it is in the same multiplet as the spin L\j2 
highest weight state Tr(Z Ll ). The operator O3 needs to be large enough to allow 
for a transition between these states which are in well separated multiplets. In 
other words, if V is too large this means that the (effect of the) operator 3 
is too small and hence the transition amplitude will be zero. A precise analysis 
would lead us to the conclusion that the three-point function should vanish for 
L' > L — N which is exactly the content of the new bound in (C.4). Note that 
Bethe states are only highest weight when the Bethe roots obey Bethe equations 
which is perfectly consistent with the requirement of on-shellness for (C.4) in 
contradistinction with (C.2). 

D Perturbation Theory for the Ising Model 

In this appendix we shall derive (4.1) and (4.2) for the partition function (3.24), 



A 




-log cx p Yl Ua f( a ~ ~ n ^ ' 



(D.l) 



{n l= 0,l} La,b=0 



with 




(D.2) 



-27- 



and ^ 

V> = ^A (D.3) 

held fixed. 

D.l Combinatoric Measure from the Partition Function 

To warm up we consider first the case where we set / — > 0. The corresponding parti- 
tion function is denoted by Z . The partition function Z just counts the number of 
permutations when the ratio of the nodes with n a = 1 and the nodes with n a = is 
fixed to 

Q= 9 —. (D.4) 
Pb 

Of course, the result is a simple binomial coefficient. In fact, this is precisely the entropy 
contribution leading to (3.12). We will now reproduce this known result in a slightly 
different way by introducing a chemical potential which ensures the correct ratio of l's 
and O's, 

Z = lim £ log / ^ I where v= P — = - 9 — . (D.5) 

a^oo A 6 ,/ 2vr x 1 p 1 + 6 v; 

The advantage of the chemical potential is that now each node becomes independent 
and we easily get 

Z = lim t log / ^ (e^ 1 -") + e *(°-*0) A . (D.6) 
a^oo A J 2n 

We just have to compute the integral over p which is saturated by the saddle point 
p = \ log 6. The above expression then simplifies to 

Z a = phgp - p a log p a - p b \ogp b = F (p) - F (p a ) - F (p b ) . (D.7) 

where 



up f p 
F (p) = p log — = / log(7r p)dp . 
e Jo 



This coincides precisely with the expected result (3.12). This is also the first term of 
the perturbative expansion anticipated in (4.2). In the next section we perform the 
perturbative expansion to higher orders. 
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D.2 Perturbation Theory 

In this section we compute the partition function (D.l) perturbatively in p. It is 
convenient to subtract from Z the combinatoric measure Z computed in the previous 
section, 

Z EE Z-Z = lim Mog^-^ t~j ~. ~ rr '- . (D.8) 

The numerator is simply given by (D.5). We will evaluate the numerator in a low 
density expansion. A small density expansion coincides with an expansion in /, see 
(D.l). For example, to leading order, the expansion of the numerator gives 



a, b " n^O 

(D.9) 

This integrals is saturated by the same saddle point H=\ log 9 as the denominator so 
that their ratio gives simply 

Repeating this calculation up to f 3 order we get 19 

^p = ^£/w + ^£/ 2 w (D.11) 

1 - Of E A") + 2 ^ E fMfMfin, -n 2 ))+ 0(f 4 ) . 
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n 

2 



3(# + l) 6 



ni,n2 



The sums can be evaluated easily at each order in p so that we can get an expansion 
in powers of p. 20 We notice that a very revealing structure arises once a derivative in 
9 is taken, 

~ 7r 2 p 3 6 - 1 7r 4 p 5 # 3 - 9 2 + - 1 7r 6 p 7 # 5 - # 4 + 9 3 - 9 2 + 9 - 1 

F)nZ ^ ——— I- - 

9 6 (# + l) 3 180 (# + l) 5 2835 (0 + 1) 7 



19 We define /(0) = to avoid writing numerous ^ 0. 
20 For example, 



6 (g+1) 2 180 + + 2835 {B+lf + U ^ > ' ^- 12 > 
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which strongly suggests that 



~ 00 f) 2n _ 1 

n=l ^ ' 

Integrating back and converting back to the p variables this translates into the following 
structure for Z (we use that for 9 = one should get Z = to integrate): 

00 

* = E ^I^ n+1 + - ^ = ^ " ^) - F( Pb ) ■ (D.14) 

ra=l 

where 

7nV 7rV 7T 6 7 

F ^ = -w--^ + ik + °^- < D15 > 

Combining this with the results of the previous section we obtain (4.2). 

By far the most important outcome of this analysis is the guess for the structure 
(D.14) of the partition function. We see that to solve this long range Ising model we 
simply need to find a single function F(p). The precise form (D.15) of the low density 
perturbative expansion of this function is not so important. In fact we derive the full 
function F(p) in section 4.1 by using a completely different approach. We find that 

F(P)= Aog Sinh(7rp W (D.16) 
Jo n P 

which matches precisely - and very nontrivially - the perturbative result derived above. 

We should emphasize that the structure (D.14) is highly nontrivial. Any small 
modification of our interaction like 

f(n) = \ log(l + p 2 /n 2 ) I log(l + p 2 /n 2 ) or \ log(l + p 2 /n 2 + 7 p 4 /n 4 ) (D.17) 

would immediately spoil the structure (D.14). This is not surprising. Integrable models 
are often isolated points in the moduli and our Ising model (D.l) seems to be no 
exception. 



E Generalizations 

In this paper we studied the quantities *4.l'(u) and B(u) in the classical limit. 21 The 
quantities Al'(vl) and B(u) are the fundamental building blocks for studying the three- 
point functions involving one non-protected operator and two protected operators. So 

21 So far we omitted the subscript L' in -4(u) denned as in (1.3) since it was always the same. In 
this section we will always keep this subscript explicit since it will change for the different examples. 
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Classical non-protected 
(9 2 = J>(u)Tr(Z...X...) 



Figure 8. A tree point function of tree classical operators. Two of the operators {0\ and 
O3) are protected and may contain SU(3) fields and fermions. The third operator, 02, is a 
non protected classical operator in the SU(2) sector. The non trivial parts of such structure 
constant are the same building blocks studied in this paper, A and B, coming from the 
normalization of 02 and the overlap between O2 and 0\. 



far, for simplicity, we considered the £77(2) setup represented in figure 2. It is simple 
to generalize this setup to more general configurations. 

Take for example the setup in figure 8. In this example the protected operators 
contain all £77(3) scalars and they can even contain fermionic operators. Still, the 
non-trivial contractions are those between operators O2 and 0\ and those are the same 
as in the previous example. Hence, we find that the ratio r defined as in (5.1) is given 
by exactly the same result (5.2), 

A v {n) 



r fig . 8 



exp 



riooo 
^123 



(simple combinatorial factor) x 



dt 



du 
2m 



q log(l — e 



iqt\ 



Luc fc 



B(u) 

— J du p log (2 sinh(7rtp)) 



(E.l) 



Another example that we can consider in the setup of figure 9. In this example we 
have an extremal correlator between one classical protected operators (O2) and two 
classical non protected operators (0% and O3). 22 There are now two non-trivial group 



2 For extremal correlators as the one in figure 9, there is also a contribution from the double trace 
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1 1 1 1 II II II 

1 1 1 1 II II II 

II II II II II 



Classical protected 
1 = £Tr(Z...X...) 



Figure 9. An extremal correlator of one classical protected operator (O2) and two non- 
protected classical operators (Oi and O3). The operators 0\ and O3 are characterized by 
the Bethe roots {vj} and {vjj} respectively. As before, the non trivial parts of such structure 
constant are the same building blocks studied in this paper, A and B, coning from the 
normalization of O2, O3 and the overlaps between these two operators with parts of 0\. 



of contractions: the contractions between O2 and 0\ and those between O2 and O3. 
Each of these contractions corresponds to the inner product between a Bethe state and 
a vacuum descendent. Let the roots of operator 0\ be {vj} and those of operator O3 be 
{wj}. Those roots will be distributed in cuts as before. We denote the cuts associated 
to the operator O a by . Similarly, the density and quasi-momenta associated to 
those cuts are, respectively, and q( a \ 23 Then 

r fig- 9 — r ,. j00 
°123 

1 

du 



C?S(v, w) = (gimple combinatorial factor) x A L v(v) x r A^(w) 



exp 

a=l,3 



dt 



2ni 



q 



(a) logfl - e iq 



B(v) B(w) 
dt I dup {a) log(2sinh(7rtp (a) )) 



UC 



(a) 



uc 



(a) 



piece of O2 that contribute at the leading planar order. Here, we only consider the single trace 
contribution. 

23 The length entering the definition (3.29) of the quasimomenta is L' for a = 3 and L — L' for a = 1. 
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Figure 10. A tree point function of tree classical SU(2) operators. Two of the operators (02 
and O3) are protected and the third operator, Oi, is non protected. For each partition of the 
roots of 0\ into left and right sets, we have the non trivial overlaps computed by A between 
O2 and 0\ and between O3 and 0\. 



Let us consider one final example, represented in figure 10. In this example, there 
is one non-protected operator 0\ which splits into two protected operators. The two 
halves of 0\ are non-trivial. This is in contradistinction to all the examples considered 
so far where all the non-protected operators were split into one trivial half and one non- 
trivial half, see e.g. figure 2. The operator 0\ is parametrized by N Bethe roots {uj} 
which parametrize the N magnon excitations on the corresponding spin chain. When 
we break that state in two, N' of these magnons will go to the right while N — N' 
magnons will go to the left. The number of roots N' is fixed by -R-charge conservation 
but we still have to sum over which magnons go to the right and which magnons go to 
the left. That is 

rag . 10 = S§M = (simple combmatonal factor) t ^-MW 

O199 _ ID ( 11 ) 

A PU0=u, |/3|=N' v ' 

To summarize, we see that the A and B are basic blocks that appear in many 
examples. Their classical limits therefore capture a general pieces of classical correla- 
tors. We like to think on that piece as the probability for a macroscopic piece of the 
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classical string to decay into a classical protected string. We call this process Classical 
Tunneling. 
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